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Abstract 

This paper presents a general framework for the design of hnear controUers for 
hnear systems subject to time-domain constraints. The design framework exploits 
sums-of-squares techniques to incorporate the time-domain constraints on closed- 
loop signals and leads to conditions in terms of linear matrix inequalities (LMIs). 
This control design framework offers, in addition to constraint satisfaction, also 
the possibility of including an optimization objective that can be used to minimize 
steady state (tracking) errors, to decrease the settling time, to reduce overshoot and 
so on. The effectiveness of the framework is shown via a numerical example. 

Keywords: Linear Systems; Constrained systems; Polynomial methods; 

LMIs; Motion systems. 

1 Introduction 

The transient response to reference commands or disturbance inputs is an important per- 
formance qualifier in many control systems. Unfortunately, most control design strategies 
cannot cope directly with requirements on time-domain signals such as actuator amplitude 
or rate limits, no output signal overshoot or undershoot, trajectory planning constraints 
and so on. Especially in the continuous-time case, there are hardly any systematic con- 
troller design methods to enforce time-domain constraints on e.g. tracking errors and 
control inputs. 

In the discrete-time case, model predictive control (MPC) (see e.g. the surveys 
is a widely used technique to cope with constraints on inputs and states. In MPC a control 
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action is prescribed that is obtained by solving a finite or infinite horizon optimization 
problem that can incorporate input, state and output constraints in a direct manner. 
A drawback of predictive control concepts and online optimization-based methods in 
general is that they require a high computational effort with the consequence that they 
cannot be implemented on fast motion systems where high sampling rates are required, 



typically in the order of several kHz. Explicit MPC |4j-[6 37 might offer an appealing 
solution as it precomputes a piecewise affine state feedback for discrete-time systems 
off-line. Still, the explicit control law often leads to a complex description consisting 
of many affine feedbacks, which also cannot realize the high sampling rates typically 
needed for motion systems of considerable size, although recent research is focussed on 
decreasing the implementation complexity of MPC, see for instance |l5|[T6||20|[23|[26] and 
the references therein. An alternative approach with strong ties to MPC is based on so- 
called reference governors, see e.g. (3||l2] and the references therein. A reference governor is 
a nonlinear device that is added to a primal controller, which functions well in the absence 
of constraints. The reference governor modifies the reference signal supplied to the primal 
controller in order to enforce the input and state constraints. This approach suffers 
from the mentioned drawbacks in MPC to some extent as well, but has the advantage 
that the reference modifications are often needed at a lower sampling frequency than the 
updates of the primal control loop. A major difference with the method presented in this 
paper is that the overall control systems in case of reference governors become nonlinear 
devices modifying the supplied reference, while the method in this paper aims at designing 
/m ear controllers that satisfy the time-domain constraints without any modification of the 
references or disturbances. 

Besides these predictive control methods, that are typically suited for a discrete-time 
context, there are only a few methods available in the literature that can directly syn- 
thesize controllers incorporating time-domain constraints in the continuous-time setting. 



For instance in the case of input constraints, 13,17 consider the linear quadratic reg- 
ulator problem with positivity constraints on the input, while various control problems 
with amplitude and rate constraints on the input signal are solved in [33j. The latter 
line of work has also been extended to stabilization and output regulation problems with 



amplitude and rate constraints on certain output variables, see, e.g., 32 . Other meth- 
ods exist that actually allow the control output to saturate such as, for instance, the 
usage of anti-windup schemes [l0|[35|[36] or LQR/LQG control methods [7, 14 . These 



methods, however, do not enforce constraint satisfaction but rather guarantee stability or 
recover performance despite the saturation nonlinearity in the loop. The above mentioned 



techniques cannot handle time-varying constraints, and, except for 32 , state or output 
constraints are not considered either. In addition, all the above mentioned techniques 
result in general in nonlinear controllers. 

As already briefiy mentioned, in this paper the objective is to derive a design method for 
linear controllers that incorporate possibly time-varying time-domain constraints on all 
closed-loop signals (inputs, states and outputs). Within this context, a commonly used 
method to capture the essence of time-domain specifications is the reformulation into 
frequency domain requirements [s]. Unfortunately, such reformulations are in general 
either approximate, conservative or both. 

A methodology to enforce time-domain constraints on the input and output of a continuous- 
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time linear control system is presented recently in where linear matrix inequality 
(LMI) techniques are used to synthesize a fixed order linear controller that satisfies the 
constraints. This is done in a polynomial setting in the sense that a controller is designed 
according to the well-known pole placement method using the Diophantine equation. 
This method allows the design of a controller that results in a closed-loop transfer func- 
tion with prescribed pole locations, either exact, or within an admissible region of the 
complex plane. In (Tol, all controllers with the prescribed pole locations are character- 



ized using the Youla-Kucera parametrization 21 . Next the degrees of freedom of the 
Youla-Kucera parametrization are used to enforce certain time-domain constraints, such 
as bounds on the input amplitude and output overshoot, exploiting sums-of-squares tech- 
niques. Unfortunately, the approach in |19j is limited to the assignment of distinct strictly 
negative real closed-loop poles, which is a severe restriction in the case of many practical 
situations such as, for instance, lightly damped systems. As a consequence, there is a 
strong need for a general framework encompassing arbitrary closed-loop pole placement. 
The development of such a framework is the main purpose of this paper. 



In particular, we propose an extension to the method in 19 , which leads to a general 
design framework based on sums-of-squares LMI techniques and we show indeed that the 
resulting linear controller satisfies the time-domain constraints on closed-loop signals, even 
when complex conjugate poles are assigned. This framework is based on two relaxations. 
One of these relaxations, of which a preliminary version was presented by the authors 
in [l], can solve the constrained control problem at hand with arbitrary accuracy and 
still lead to LMIs. In addition to constraint satisfaction, we will also include an objective 
function in the convex programming problem that can be used to minimize steady state 
(tracking) errors, to decrease the settling time, to reduce the overshoot and so on. As a 
consequence, the ideas presented in this paper will lead to a general design framework for 
optimized linear controllers with guarantees regarding constraint satisfaction. 



The organization of the paper is as follows. The proposed methodology from 19 is briefiy 
reviewed in Section [2j The extension to complex conjugate poles is treated in Section [3| 
which includes the main results. Section |4] discusses the proposed control design method, 
and Section [5] provides an illustrative example. Finally, the conclusions are stated in 
Section [6l 



2 Methodology involving real poles 



In 19 a method is presented to incorporate time-domain constraints on input and output 
signals of a linear system. It is shown that finding a controller of fixed order that satisfies 
these constraints boils down to solving a set of LMIs. In this section, we shortly review 
this procedure for completeness and self-containedness. 
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2.1 Youla-Kucera parametrization 

Consider the control system depicted in Fig. [T] with a hnear single-input-single-output 
plant P given by the strictly proper transfer function 



Pis) 



b{s) 



ais 



(1) 



where a{s) and b{s) are polynomials in the Laplace variable s. The controller C, which 
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Figure 1: Block diagram of the closed-loop system with controller C, plant P, and refer- 
ence signal r, control output signal u, and output signal y. 



is to be designed, is described accordingly by 



C(s) = % (2) 

c(s) 

resulting in the complementary sensitivity given by 

T( ) = ^ = Ks)d{s) 
^' r{s) a{s)c{s) + b{s)d{s)' ^' 

If a{s) and b{s) are coprime (i.e., their greatest common divisor is 1), then arbitrary pole 
placement can be achieved by designing the corresponding controller polynomials. This 
is done by solving the polynomial Diophantine equation 

a{s)c{s) + b{s)d{s) = z{s), (4) 

where z{s) = (s +Pi)(s +P2) ■ ■ ■ {s +p„) is the polynomial with given roots —pi, . . . , —Pn, 
which are the desired poles of the closed-loop system. There are infinitely many solutions 
(c(s),(i(s)) to Q, but there is a unique solution pair (co{s) , do{s)) such that deg do(s) < 
dega(s). In this case we have that do{s) is of minimal degree and as such, {co{s) , do{s)) 
is called the ci-minimal solution pair. All possible solutions to the Diophantine equation 
can then be written as 

c{s) =Co{s) + b{s)q{s), 
d{s) = do{s) - a{s)q{s), 

where q{s) is an arbitrary polynomial such that co(s) -|- b{s)q{s) is non-zero. This poly- 
nomial, called the Youla-Kucera parameter |9], creates extra freedom in the design of 
the controller. While the closed-loop poles are invariant for any choice of the Youla- 
Kucera parameter, the Youla-Kucera parameter enables placement of closed-loop zeros to 
alter the response. Only proper controllers are considered and therefore there is a degree 
constraint on q{s). Since the plant was assumed to be strictly proper, and under the 



4 



additional assumption that deg2;(s) >2dega(s) — 1 (to enable arbitrary pole placement 
with proper controllers), this constraint is given as in |22j by 

degg(s) < deg z{s) — 2dega(s). (6) 

The extra freedom in the control design parameterized by q{s) satisfying ^ can now be 
used to satisfy additional time-domain constraints as will be explained in the next section. 



2.2 A positive polynomial formulation of time-domain constraints 



We will explain the procedure in [19] using the typical example of constraints on the step 
response. Hence, we consider the response y to a step input (r(s) = ^). The Laplace 
transform of the closed-loop system's output (assuming zero initial conditions) is then 
given by 

lb{s)d{s) lb{s)do{s) la{s)b{s) 

y{s) = = -—q{s). 7 

S Z{S) s z[s) s z{s) 

At this point of the control design a restrictive assumption was made [l9], namely 

Assumption 2.1 All the assigned poles —pi, . . . , — p„ are distinct strictly negative ratio- 
nal numbers. 

Using this assumption and z{s) = Y[^=iis+Pi) the partial fractional decomposition of ([T]) 
leads to 

vM = ±^, (8) 

where po = and yi{q), i = l, . . . ,n are appropriate coefficients following from the decom- 
position, which are influenced by the choice of the design parameter q{s) = Yltlo The 
coefficients yi{q) depend in an affine manner on the parameter q = (go, qi, ■ ■ ■ , qdq) in the 
sense that there exist matrices A e ]R(«+i)x("+i), B e M("+i)><('?d) and a vector b e M"+^ 
such that 

/ yo{q) \ 



A 



Bq' + b. (9) 



\yn+i{q)' 

This follows directly by comparing ([T]) and (|8|, and equating the coefficients of the powers 
of s in the resulting numerator polynomials (see also (46) below for an example). The 
corresponding time-domain signal is given by 

n 

y{t,q) = Y.y^{q)e-^''- (10) 

1=0 

Let Pi = ^ be the ratios of the integers Ui and di, and let m denote the smallest positive 
numbeinof the denominators such that Pi = ^ for some positive integers pi, i = 0,1, ... ,n. 



^In principle m can be chosen to be any positive number that resuhs in integer values of pi. However, 
it turns out that by choosing m as the smallest possible positive number the order of the resulting 
polynomial optimization problem is the lowest. 
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This means that the time-domain output signal at time t G M+ := [0, oo) can now be 
expressed as the polynomial 



y 



(A,g) = X^l/.(g)A« (11^ 



i=0 



in the indeterminate A = e Obviously, A lies in the interval [0, 1] as tGM+. Suppose 
that the output y(t, q) of the system needs to be bounded according to 

l/min < y{t, q) < 1/max V t G IR+. (12) 



Formulation (|12|) is equivalent to enforcing the polynomial bound constraints 

P,(g,A):=,(A,g)-, |0 vAg[0,1], (13) 
^2(g, A) := ?/max - 2/(A, g) > 

where Pi and P2 are polynomials in both A and q. This problem is a special case of the 
following more general problem of minimizing a polynomial with polynomial constraints 
over a basic semialgebraic set. 

Definition 2.1 A set D is called a basic semialgebraic set if it can be described as 

V = {xeW \ ei{x)>0, i = l,...,Me and 

/,(x)=0, j = l,...,M;} ^^^^ 

for certain polynomials Ci : M"^ — J-M, i = l, . . . , and fj : M"^ — >]R, j = 1, . . . , Mf. 

Problem 2.1 (Polynomial optimization problem) Consider two variables z G M"^ 
and X G M""^ and let polynomials gt : M"^ x M"^ — )■ M, i = 1, . . . , M^, and p : M"^ — M 
6e given. Moreover, let a collection of basic semialgebraic sets Vi C M"^, / = 0, . . . , be 
given. A (robust) polynomial optimization problem according to this data is given by 

min p{z) 

N (15) 
s.t. gi{z,x)>0, i = l,...,Mg y xe[jVi. 

1=0 



Indeed, (13) can now be written in the form of Problem 2.1 by taking z = q, x = X, M, 



'9 ■ 



2, N = 0, piz)=0, gi{z,x) = Piiq,X), g2{z,x) = P2{qA), and I^o = {A G M | < A < 1}. 



Although the bounds ?/min and i/max in (13) are chosen to be constants for illustrating 



purposes, they can also be selected as polynomials in A, i.e., in the form ?/min(A) and 



l/max(A) without any complications. In this case the bounds in (12) become time-varying. 



Univariate positive polynomial constraints (meaning polynomials in only one variable). 



such as (13) with A G X'o = [0, 1] ^ M, can be transformed into LMI conditions, see [19 



for the details. Once we transformed the design problem into a polynomial optimization 



problem as formulated in Problem 2.1, there are appropriate tools available for solving 



the problem. Therefore, we restrict ourselves to the transformation of the constrained 



control problems at hand into manifestations of Problem 2.1 
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The approach discussed in this section is not restricted to bounding only the output of a 
system. Indeed, by using the appropriate transfer functions, any signal in the loop can 
be constrained. The control output u, for example, can be bounded using 

, , , .a(s)do(s) , ,a^(s) , , , , 

= ^(,) _ r{s)^q[s) (16) 

in addition to, or instead of ([T]). Also, other Laplace transformable inputs or disturbances 
can be used as long as the poles of the Laplace transform of the corresponding signal are 
distinct strictly negative rational numbers and differ from the closed-loop poles p^. In case 
the disturbance signal is a (filtered) random process, the method cannot be applied as 
is. However, a possible extension can be to bound the infinity- or 2-norm of the suitably 
weighted process sensitivity (or other relevant transfer functions) via e.g. 



This way the knowledge of stochastic disturbances can be used to shape the relevant sen- 
sitivity functions to achieve desired disturbance reduction. 

A combination of requirements on different reference signals can easily be handled at the 
cost of increasing the size of the set of LMIs. As for LMIs there are efficient solvers 



available, e.g. 34 , transforming the problem at hand into Problem 2.1 provides an ef- 
fective solution. The problem derived in this section was only a feasibility problem (the 
cost criterion -piz) was chosen to be 0). In Section [i] below, we will also provide relevant 
choices for the cost criterion that next to satisfaction of the time-domain constraints also 
provides additional desirable properties of the constructed controller. 



3 Problem formulation: the complex poles case 



The polynomial representation (10), as derived in 19 , of the time response of a linear 



system to a Laplace transformable input is unfortunately only possible when strictly 



negative rational closed- loop poles are assigned (see Assumption 2.1). However, in many 
cases the assignment of purely real poles can be undesirable, especially in lightly damped 
systems such as most motion systems. Furthermore, many reference signals have Laplace 
transforms with complex poles. If, for instance, a sinusoid is used as the reference signal 
instead of a step, the Laplace transform is given by r(s) = ^5^^ resulting in complex 
poles in the system's response. Therefore, such reference signals cannot be handled by 
the approach from The main objective of this paper is to present a solution to the 
linear control design problem for time-domain constrained systems of which the Laplace 
transforms of the closed-loop responses may contain complex roots. 
When we allow both distinct real and complex poles to be present in the closed-loop 
transfer function T(s) and/or the Laplace transform of the reference signal r(s), the 
Laplace transform of the system's output can be decomposed as the partial fractional 
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decomposition 



i/(^)=E — 



nr+nc/2+l 

+ y '-^ + 

Z-^ Q 4- -I- i R- 



i=nr+l ^ + + ^/^i s + ai- j(3i 

where rij. and Uc denote the number of real and complex poles, respectively, —pi, i = 
1, . . . ,nr are the locations of the real poles, — i = rir + nc/2 + 1 are the locations 
of complex conjugate pairs of poles, and j/j are the possibly complex coefficients (with 
complex conjugate y*) that affinely depend on the design parameter q, see ^ (we omitted 
this dependence on q for ease of exposition). To enforce stability, we again assume that 
the assigned closed-loop poles have strictly negative real part. The corresponding time- 
domain signal is then described by 

y{t) = Y,yre^'"' + E (y.e-^'^'* + l/*e^'^'*)e-"'*. (19) 

As before, we use the following assumption 
Assumption 3.1 pi, ai, and Pi are rational numbers. 

We denote p,- = — , a,- = — , r = — where m is the smallest positive number (not necessarily 

an integer) of Pi and such that Pi and cij can be taken as integers. We denote /3j = — 
for a number Q (not necessarily an integer) such that /3i can be taken as integer as well. 



For guidelines how to choose ^, see Remark 3.1 below. Furthermore, let 



A = e-\ (20) 

Using Euler's formula e^'^ = cos(0) + j sin((/)) and decomposing the complex coefficients 
as yi = ai+jbi, y* = ai-jbi, yields 

rir 

y{t)=Y,y,y^ 

1=0 

nr+nc/2+l y ) 

+ E (^^2 cos(^i^r) + bi2 sin(^,^r)) A"" . 

i=nr+l 

Obviously, the terms involving the complex poles are non-polynomial in the indetermi- 
nate A because of the presence of cos(/3j6'r) and sin(/3j^r), which make it impossible to 
directly use the positive polynomial approach in Section [2] to bound the output as in 
( [l2| ). Although the parameters (3^ and pj are fixed as a result of the pole placement 
(and the choice of m and 6), there still is freedom in the choice for the coefficients aj,6j, 
which depend on the coefficients q = (go, • • • ,5^,) in the Youla-Kucera parameter q{s). 
We propose two relaxations to determine the values a,, bi via polynomial optimization 



problems to shape the time response y{t), thereby overcoming the limitations in 19 . 
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The first approach is based on an exponential bound relaxation that results in univari- 
ate polynomials. This method has the advantage that it results in a simple polynomial 
optimization problem, but introduces some conservatism. The second method proposes a 



multivariate polynomial relaxation that leads to polynomial problems as in Problem 2.1 



while the conservatism can be made arbitrarily small. Both methods result in polynomial 



optimization problems of the type as in Problem 2.1 that can be solved using LMIs 



3.1 Exponential bounds relaxation 

To resolve the problem induced by the presence of the products cos(/3j6'r)A"' and sin(/3j6'r)A"' 



in (21) we relax the problem by using the fact that 

cos(^i^r), sin(^i^r) G [-1,1] Vr G M, (22) 



and instead of the exact time- response (21), we consider 



rir nr+nc/2+1 

l/upper(A) = Y^ViX^^ + J2 (2|a,| + 2\k\) X"', 

1=0 i=nr+l 
Hr nr+nc/2+l 

Z/iowcr(A) = 5^i/.A^' - (2|a,| + 2\h\) A^\ 

i=0 i=nr+l 



(23) 



In contrast to (21), these exponential bounds on the closed-loop time response are uni- 
variate polynomials in the indeterminate A = e~'^ (if q is fixed) and can be bounded by 
specified polynomials 5'u(A) and gi{\) via the polynomial non-negativity constraints 

^3(g, A) := ^„(A) - |/uppcr(A) > ^ ,n .^.s 

P4(g,A):=yiower(A)-^7KA) >0 ^^^^^'^J' ^^^^ 

where we included the explicit dependence of and hi on q again. The constraints 



(24) cannot straightforwardly be cast in the form Problem 2.1 because of the nonlinear 
operator | • |, which is present in these equations. However, each of the two nonlinear in- 
equality constraints P3 and P4 can be expressed as 2^"^^ equivalent polynomial inequality 
constraints P3 and P4 (2 inequalities for each absolute value expression). Enforcing non- 



negativity of (24) on the interval A G [0, 1] is then again a special case of Problem 2.1| with 
z = q, x = X, Mg = 2, Mh = 0, N = 0, p{z,x) = 0, gi{z,x) = P-i{q,\), g2{z,x) = Pi{q,\), and 
"Do = {A gM I < A < 1}. Therefore, it is possible to determine the values q = (go, • • • , Qdq) 



such that the upper and lower bounds (23) of the closed- loop time response are bounded 
by guiX) (e.g. 5'„(A) = y„,ax) and gi{X) (e.g. guiX) = y-min) via a polynomial optimization 
problem. The exponential bounds relaxation does introduce some conservatism by using 



relaxation (23) instead of the exact time- response (21). The second method presented 
next offers the possibility to render this conservatism arbitrary small. In other words, 
the second method can approximate the original time-domain constraints with arbitrary 
accuracy and still lead to polynomial optimization problems. 
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3.2 Multivariate polynomial relaxation 



The time response (21) is equivalent to 

nr+nc/2+l 



J2 (cos(A^r) -J sin(/3,0r)) 



(25) 



i=nr+l 



-Ha^-jh) (cos(/3,^r)+jsin(/3,^r))] A°'. 

De Moivre's formula, which is closely related to Euler's formula and {e^'f')'^ = e^"""^, states 
that for any G M and any integer n G Z 

(cos(0) + j sin(0))" = cos{n(f)) + j sm^ncj)), (26) 



and hence (25) is equal to 

rir 



i=0 
nr+nc/2+1 



(K + J^^) [cos(0r) - jsin(er)]''- + 



(27) 



{ai - jbi) [cos(^^r) + J smi9T)f^ 



Appropriate polynomial functions Wi : 



R and r, : _^ 



Ur + I, . . . ,11^. + 



nc/2 + 1 in two variables can now be defined such that (27), and thus the time response 
(21), can be written as 



(28) 



1=0 

nr+nc/2+1 

(aj2wj(cos(6'r), sin(6'r)) 

i=nr+l 

+6i2ri(cos(er), sin(^r))) 
This proves the following theorem. 

Theorem 3.1 Consider the closed-loop system ^ and let y be the response to a reference 
input r and assume that the Laplace transform y{s) of y has only distinct poles such that 
(18) and Assumption 3.1 hold. Then we have that 

{y{t) I t G M^} = {y{u, V, A) | (m, V, A) G J^original}, 

where y{u, v, A) is given by the multivariate polynomial 

i=0 

nr+nc/2+1 

+ {ai2wi{u,v) + bi2ri{u,v))X''' 

i=nr+l 



(29) 



(30) 
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with : — 7- M and : — ?■ M, z = ra^ + 1, . . . , + nc/2 + 1 polynomials as in (28 ) 
and 

^original ■= {{u,V,\) E \ U = COs{9t) , V = Sm^Or) , 

X = e ^ for some r G M | . 



Proof: The reasoning before the formulation of the theorem revealed that y{t) under 



the given assumptions is equal to (19), which can equivalently be written as (28), where 
A = e"'^ and r — * ~ 



m 



Due to Theorem 3. 1 , bounding the output as in ( 12 ) to the interval [ymin, Umax] is equivalent 
to enforcing the polynomial non-negativity constraints 

Pbiq, u, V, A) := y{u, V, A) - ymm > 0, .^2) 
Pe{q,u,v,X) := y^^^-y{u,v,X) >0 

for all (u, f , A) G J\)riginai- Recall that y{u,v,X) depends on q via yi,ai and 6,. As we 
mentioned before, it is of interest to transform the linear constrained control problem 



into Problem 2.1 The conditions (32) are not in this form due to the fact that J\3rigmai 



is not a (finite union of) basic semialgebraic set(s) as in Definition 2.1 However, this set 
can be overapproximated by a finite union of basic semialgebraic sets in an arbitrarily 
close manner. 

Definition 3.1 We call a set J^approx on e-close over approximation of original for some 
e > 0, if it satisfies the following three properties: 

1- approx = Uzio f^'^ ^ finite collection of basic semialgebraic sets J-q, . . . , J-jv; 

2. J~ original ^ approxi 

3- J^approx C J^original + ^e, whcrC := {(0, 0,z) \ S < Z < E. 

Hence, an e-close overapproximation of J-'originai contains the set J-'originai as drawn by the 
white line in Figure |2] (for 9 = 1), but it is e-close in the sense of property 3. Hence, for 
small e > 0, replacing J-'originai by J-kpprox only results in small errors and all guarantees on 
-^approx also apply to J-'originai due to property 2. Moreover, due to property 1 an e-close 
overapproximation v^approx of J-'originai can be used to embed the polynomial constraints 



in (32) for all (u, f , A) G J^ariginai into a version of the constraints in Problem 2.1, where 



Vi = J^i,l = 0,...,N. 

The following algorithm provides an algorithm that constructs for each desirable level of 
approximation e an e-close overapproximation of J^^riginai- The basic idea of the algorithm 
is to overapproximate the J^origmai-set by the union of basic semialgebraic sets J^i, which 
are obtained by sphtting the set J-'originai in the r-direction by considering intervals Ii : = 
[tuti+i), I = 0,...,N, where = tq < ri < . . . < tn+i = oo. On each of these 
subintervals // we approximate e""^ by '?/^i(cos(^r), sin(^^r)) using Fourier series, where 
: — M is a polynomial such that je""^ — iI!i{cos{6t) , sin(^r))| < e for all t E Ii. Next 
to e, the algorithm uses another parameter 0<T< which indicates the desired length 
of the intervals / = 0, . . . , — 1 (although it will be modified such that all intervals 
have the same length). 
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Figure 2: J-'originai (white line) drawn inside the cyhnder given by u'^+v^ = 1 and < A < 1. 

Algorithm 1 Let 0<e<l andO<T<^ be given. 

Step 1: Define N -.— l^^l and tn ■— — in. e and tn+i ■— oo . 

Tn := {{u, v,X) eR^ \u'^ + ^1 andO<X< e}. (33) 

Step 2: Divide the remaining interval [0,tn) in N subintervals of length T:—^ <T < 
f. Ii:^[Ti,Ti+i) withTi^lf,l^O,...,N-l. 

Step 3: For each l — 0,...,N — l define a function 0^ : R — >■ R that satisfies: 

• (pi is at least continuously differentiable, but preferably m times continuously 
differentiable (C"^) for m & N large; 

• (pi is periodic with period 

• (/)i{r) = e^^ for all r E Ii. 

Step 4: For each I = 0, . . . , N — 1 compute the Fourier series approximation of (pi of 
sufficiently high degree Ki such that 

Ki 

\<Pi{t) - ^[afc cos{k9T) + 6fcSin(A;6'r)]| < e for all r e //, (34) 

fc=0 

where ak,bk, k = 0, . . . , Ki are the Fourier coefficients of (pi. 

Step 5: For each 1 = 0, ... , N—1 use De Moivre's formula to rewrite X^^qI^*; cos(A;^t) + 
bk sm{k9T)] obtained in the previous step as 

Ki k 

5]^Cfc,(cos(^T))*^(sin(^T))' =: V'Kcos(^t), sin(^T)), 

A;=0 i=0 

where : R^ — >■ R is a polynomial of degree equal to the degree of the Fourier series. 
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Step 6: For each 1 = 0, N —1, define 



J^i ■= {(u, v,X) eR^ \ -e < X- ^i{u, v) < sA 

5;+iQ-Q+i5;<0}, (35) 

where Ci := cos{6ri), Si := sva.{6Ti). 
Step 7: Take J- approx — U/=o'^'' 

Theorem 3.2 For each < e < 1 and < T < ^ Algorithm i] produces an e-close 
overapproximation T approx of T original in the sense of Definition 3.1. 



Proof: First of all, we write J'originai as IJ^g -?originai,« with ^original,; := { (cos 6'r, sin 6'r, e"^) 
r G /;} for / = 0, 1,...,A^ as Uilo-^' ~ [0,oo). Step 1 considers the interval 1^ '■ = 
[tn,tn+i) = [—In £,00) for which it holds that < e"'^ < e. Hence, clearly J-'origmai,Af ^ 
J^N and J-Tv ^ ^originai,Ar + '^e- The construction of functions (pi in step 3 is possible 
as e""^ is continuous and the fact that < T < T < Hence, step 3 can always be 
taken, while the function 0; can still be made ^-periodic and continuously differentiable. 
Step 4 can be realized, because the Fourier series converges uniformly to a continuously 
differentiable periodic function, see, e.g., |2,29 . Therefore, uniform convergence proves 



the existence of a finite Ki such that (34) holds. Note now that for {u, v, A) G J^^rigmai,; it 
holds that u^ + v"^ = 1. Obviously, (cos(6'r), sin(6'r)) for r &Ii = [r^, r^+i) lies in one of the 
half spaces generated by the straight line in through the points (cos(6'r/), sin(^ri)) and 
(cos(6'r,+i),sin(6'r,+i) given by 

[sin(6'r/) — sin(6'ri+i)] cos(6'r) +[cos(^^ri+i) — cos(6'r/)] sin(6'r) 

+ sin(^ri+i) cos{9ti) - cos(^r/+i) sin(^^ri) = 0, (36) 
see Figure [3] In particular, for all {u, v, A) G J^originai,/ it holds that 

{Si - Si+i)u + (a+i - Ci)v + Si+iCi - Ci+iSi < (37) 



where Ci :=cos(6'r;). Si ■.= sm{6Ti). Due to (34), step 5, and using the above observations, it 



holds that J-'origmai,; '^J^i- Moreover, similar reasoning using (34) shows that J-'i'^J-'origma.i,i+ 
Ms- Hence, by taking J-'approx as the union of the resulting sets, i.e. ^approx = Ujlo -^^5 
e-close overapproximation J-'approx of J^originai is obtained. This completes the proof. ■ 

For instance, for 6' = l,£: = e~^'^'^~0.009 and T = T = 0.757r the e-close overapproximation 
J~ approx of J~ 3rig inai can be generated with tq = 0, Ti = 0.757r, T2 = l.Svr and = 00 and 
polynomials 



^o(M,t;)=0.398M-0.971f +0.616M^-0.192Mt; + 1.179t;^ 

(3^ 

+ 0.013Mt;^-0.011t;''. 



-0.015^2+0.184^2^;, 
tP^lu, i;) =0.033u+0.096t;+0.0760M2 + 0.0534Mt; + 0.094t;2 
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cos{8r) 



Figure 3: Top- view of the cylinder {{u, v, A) G | u^+v"^ = 1, < A < 1} (solid grey), the 
hyperplane {{u,v,\) G M.^ \ (36)} (dashed), and the resulting part of the cylinder given 



by {{u,v,X) G I u^+v^ = l,l^,0 < A < 1} (solid black). 
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Figure 4: Functions e ^ (solid black), ipQ (solid grey), and ipi (dashed grey). 



This overapproximation and the polynomials are illustrated in Figure |4j If one is sat- 
isfied with an overapproximation accuracy of e = e~^'^^ ~ 0.009, then one can use this 
precomputed overapproximation (for the case 6 = 1). If it is desired to have a simpler 
overapproximation (with less regions and polynomials of lower degree) or an even tighter 
approximation with e < 0.009, one can run Algorithm [T] to obtain it. 

Remark 3.1 A few comments are in order: 

• The reason to take the functions (pi in step 3 as smooth as possible (m in as large 



as possible) is that a lower degree Ki is needed to satisfy (34). Indeed, when cpi G C 



then the Fourier coefficients satisfy k^at — ?■ 0, k^hk when A; — > oo and thus the 
approximation error \ X^fc^ot'^fc cos(fc6'r) +6^ sin(fc^r)] — e^^| = | Y1^=Ki t'^fc cos^kOr) + 
bk sm{k6T)]\ on Ii is smaller than e for smaller values of Ki. 



A sufficiently high degree of Ki such that (34) holds can be obtained by increasing 



incrementally until (34) is satisfied. If one is satisfied with an overapproximation 



accuracy of e = e , then the precomputed overapproximation (38) with N = 2 
and Ki = 3 can be used in case 6 = 1. 



14 



Instead of selecting T a priori, we can also select a maximal degree K of the approx- 
imation functions ipi in the sense that Ki < K for all I = 0,1, N — 1. Instead 
of increasing the degrees of the approximation functions ipi, one now can split the 
time interval [0, — Ine) into smaller pieces until lipiicosOr, sinOr) — e^'^\ < e for all 
T E Ii is satisfied for the fixed (low) degree K . This might lead to more regions 
(a larger N). This indicates that there is a trade-off between N (number of basic 
semialgebraic sets in the overapproximation J^approx) O'^d K (the maximal degree of 
the Fourier series approximation). The smaller N the higher K and vice versa. 
However, note that the example of the overapproximation with N = 2 and K = 3 
already provides a very tight approximation of e = e~^'^'^ ^0.009 in case 6 = 1. 

The reason for splitting the set J^originai the t = ^ direction is that the exponential 
function A = e~'^ can generally not be e-close approximated with the basis functions 
cos(^r) and sin(6'r) in the interval [0, — Ine] (unless e is chosen relatively large or 9 
is very small and thus T can be selected larger). The additional scaling 9 can be used 
to adjust the period of cos{9t) and sin(^r) thereby offering a trade-off between the 
required amount of intervals N and the order of the polynomial y{u,v, X ) (p rovided 
the integer requirement on /3j is fulfilled). Indeed, N = 1 in Definition 3.1 can be 
obtained by choosing 9 sufficiently small such that e"'^ can be e-close approximated 
using u = cos{9t) and f = sin(6'r) in the interval [0, — Ine] for any e. However, this 



implies that the polynomial order (3^ (see [27]) increases due to the relation /3, 



which is undesirable from a computational point of view. On the other hand, if 9 
is chosen as large as possible (while respecting the integer requirement on /3^), the 



order of the polynomial y{u,v. A) in (27) is minimal but a large amount of regions 
N could be required. 



Using the e-close overapproximation J-'approx, the conditions (32) again are a special case 
of Problem |2T with z = q, x = {u,v, X), i = 2, m = 0, p{z, x) = 0, gi{z, x) = P5{q, u, v. A), 
g2iz, x) = PQ{q, u, V, A), and T>i = Fi, 1 = 0, . . . , N. However, since the polynomial positivity 
constraints ( p2) ) are now multivariate (meaning polynomials in more than one variable) 
instead of univariate, the equivalent LMI expression from 19 is not applicable. In gen- 
eral, checking positivity of a multivariate polynomial on a basic semialgebraic set is a hard 
problem, but it can often be approximated as closely as desired by a hierarchy of con- 
vex relaxations 24,25 . In this case, to make Problem 2.1 computationally tractable, the 



inequality conditions from (15) are replaced by stronger conditions in terms of primal mo- 
ment and dual sums-of-squares (SOS) problems to formulate a hierarchy of upper bounds 
on the minimum in Problem 2.1 that converge in the limit to the real minimum. We refer 
the reader to 24 , 25 for the conversion techniques to obtain the LMIs and further details. 

, SOSTOOLS I30i or YALMIP |27l, 



There exist software packages, such as GloptiPoly I 
that automatically build up a hierarchy of LMI relaxations, whose associated monotone 
sequence of optimal values converges to the global optimum. Numerical certificates of op- 



timality are also available, in terms of ranks of embedded moment matrices, see 24,25 
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4 Control design procedure 



To summarize the previous design setup, suppose that the c?-minimal controller (|2| for 
system ([T]) has been designed via the Diophantine equation ^ such that the desired 
closed-loop pole locations are achieved. Also suppose that more closed-loop poles are 
assigned than twice the number of poles of the plant. Then, according to (|6| there is ad- 
ditional control design freedom parameterized in the form of the Youla-Kucera parameter, 
which can be used to shape the closed-loop time response. Time-domain constraints can 
be imposed using one of the proposed relaxations in Section |3| which yield the polynomial 



constraints in (15). Next to constraint satisfaction, (15) allows also the minimization of 
an objective function p{z). This minimization can be exploited to obtain additional de- 
sired properties of the response in terms of the design parameters q and thus yi, ai and 



hi. Consider the step response for example, which can be written as (18) with po = and 
where ?/o is the steady-state solution. Desirable properties of the unit step response are, 
for instance, a zero steady-state error, a small settling-time and small overshoot. These 
properties can be accommodated in p{z) as follows 

Small steady-state error: Set p{z) = {l — yQY to minimize the steady-state error. 
Short settling-time: Set p{z) = af + bf, where index i corresponds to a slow mode in 



(21), to minimize the contribution of this mode, which improves the settling time. 
Alternatively, exponentially decreasing constraints can be specified that directly 
impose a certain desired settling behavior. 

Overshoot minimization: In case one is interested in constraining the response by a 
fixed constant, e.g. constraining the overshoot, the exponential bounds relaxation 
is not suitable. This is due to the fact that the peak values of the systems response 
often occur in the time interval where the exponential upper and lower bounds are 
still very far from the actual signal. As opposed to the exponential bounds relax- 
ation, the multivariate polynomial relaxation from Section 3.2| is typically suited 



to minimize the overshoot of a step response by constraining the response (30) as 
y{u, f , A) <7 and specify = 7 to minimize the overshoot. 

Of course, one can combine the above objectives in p{z) using suitable weighting factors. 
Furthermore, the proposed relaxations enable the incorporation of the extensions that 



were mentioned in Section 2.2 (e.g. related to responses to disturbances) in case the 



poles of the Laplace transform of the corresponding signal are complex. Setting up the 



optimization problem ( 15 ) by including the time-domain constraints on closed-loop signals 
using one of the proposed relaxations and defining the objective function p{z) provides a 
systematic manner for obtaining linear controllers with desirable properties. This design 
framework will be illustrated in the next section. 
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5 Numerical example 



We start with a simple simulation example to illustrate the efficiency of the proposed 
design method. Consider the simple model given by 

Pis) = 44 = (39) 
^ ^ u{s) s + 1 ^ ^ 

The control objective is to let y track a step reference from to 1 as close as possible. 
Moreover, the controller ^ will be designed such that the assigned complex closed-loop 
poles are pi^2 = — 1 ± 2j, pa 4 = — 2 ± 4j. This is done by solving the Diophantine equation 
Q leading to the (i-minimal controller 

_ rfo(s) _ 68 

cois) s3 + 5s2 + 28s + 32' ^ ' 

resulting in the closed-loop system given by the complementary sensitivity function 

T(s) = 5 . (41) 

^ ' + 6s3 + 33s2 + 60s + 100 ^ ' 



Using the Youla-Kucera parameter q(s) and realizing that according to Q we have 
degg(s) < 2, i.e., g(s) = qo + Qis + q2s'^ , the set of allowable controllers assigning the 
specified closed-loop poles is parameterized as 

^^^^ d{s) do{s) - a{s)q{s) 



c{s) coi-s) + b{s)q{s) 
^ 68 - (go + (go + qi)s + (gi + ga)^^ + 92^^) 
s3 + 5s2 + 28s + 32 + (go + gis + g2s2) ' 

resulting in the set of closed-loop transfer functions 

X ^ 68 - (go + (go + gi)s + (gi + q2)s^ + gas^) 

s4 + 6s3 + 33s2 + 60s + 100 ' ^ ' 



The Laplace transform of the time response of (43) to a step input is then parameterized 
as 



vis) = -Tis) = 

68 - (go + (go + gi)g + (gi + ga)^^ + q2s^) 
s(s + 1 + 2j){s + 1 - 2j)(s + 2 + 4i)(s + 2 - 4j) ' 

The corresponding partial fractional decomposition is equal to 
y{s) =— + , , , ^. + —— — + 



(44) 



s s + 1 + 2j s + 1 - 2j s + 2 + 4j 
s + 2-Aj 



(45) 
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where yo, ai, 61, 02, &2 can be solved from the hnear system of equations 
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(46) 



where qo,(li,(l2 are the free variables in the Youla-Kucera parameter to shape the time 
response. The goal is to determine values of yo,ai,bi,a2,b2 (via qo, (11,(12) such that the 
closed-loop time response to the step input has a favorable shape. The LMI problems 
were modeled with YALMIP 



27 and solved with SeDuMi 34 



5.1 Using the exponential bounds relaxation 



The exponential bounds on the step response of the closed-loop system (43) are given by 

^uppcr(t) =yo + (2|ai| + 2|6i|)e-* + (2|a2| + 2\b2\)e-^\ 
yiowcr(t) = 2/0 - (2|ai| + 2|6i|)e-* - (2|a2| + 2|62|)e-'*, 



where j/o, ai, 61, 02, &2 are related to qo,qi,q2 via (46). The goal of this relaxation is to 
determine go, ^i, q2 such that 

= 9u{X) - 2/uppcr(A) > 0, , . 

PaW = mo^orW - gi{X) >0, ^^^^ 

for appropriately chosen time-varying hounds related to (7u(A) and gi{X). Note that A = e~* 
and 



2/uppcr(A) =yo + (2|ai| + 2|6i|)A + (2|a2| + 2|62|)A' 
2/iower(A) =yo- (2|ai| + 2|6i|)A - (2|a2| + 2\b2\)\^. 



(49) 



To define the time-varying bounds, we first consider suitable upper and lower bounds on 



the step response that correspond to the ci- minimal controller (40) (with g = 0). These 
are given by 

gu .. ,(A) = 0.68 + 1.58A + 0.38A^ 

"original V / ' ' " f Vir\\ 

c/i . . ,(A) = 0.68 - 1.58A-0.38A^ ^ ^ 

'■original \ / ' 

respectively. Based on these bounds, we now define tighter bounds guW and gi{X) that 
are specified to be 



9u{\) = ^7«,.,„.,(A) + c„ = 1.01 + 1.58A + 0.38A' 
^?/(A) = ^?Wina.(A) + Q = 0.99 - 1.58A - 0.38A 



(51) 



where 0^ = 0.33 and c/ = 0.31 to guarantee a small steady-state error (smaller than 0.01). 



The dominant term in the upper and lower bound in (49) corresponds to the slow mode 



e * with coefficient (2|ai| + 2|6i|). To minimize the contribution of this term and to 
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bring the steady-state error close to zero, we specify the objective function according to 
Section |4] as 

p{qo, gi, gs) = 10(1 - yof + 2(a? + bj), (52) 



where yo, ai and bi depend on go, gi, 12 as in (46 ). This results in the optimization problem 



min p(go,gi,g2) 
go ,91, '72 

s.t. 



(46) 



(53) 



(48) VAg[0,1]. 



The Youla-Kucera parameter resulting from the minimization problem (53) is 

2 



q[s) 



-32.0 - 23.0s - 3.0s' 



which yields the controller and closed-loop 

3s^ + 26s2 + 55s + 100 



Cis) 



T(s) 

together with the new bounds 



s3 + 2s2 + 5s 

3s3 + 26s2 + 55s + 100 
s4 + 6s3 + 33s2 + 60s + 100' 



^7n_(t) = 1.00 + 1.25e- 
^7^„ew(i) = 1-00 -1.25 



-2t 



-2t 



(54) 

(55) 
(56) 

(57) 



This shows that the steady-state error is zero and the contribution of the slow mode 
has been completely eliminated. The step responses together with their bounds of the 



original closed-loop system (41) and of the new, optimized closed-loop system (56) are 
depicted in Fig. [sj^a), while the Bode diagrams of the original controller (40) and the new 
controller (55) are shown in Fig. [sj^b). From Fig. |5](a) it is obvious that the step response 



of the designed closed-loop system satisfies the specified bounds and additionally results 
in zero steady-state tracking error. Upon examination of Fig. Islb) this can be explained 



by the fact that controller (55) exhibits an overall higher gain and hence results in a 



higher bandwidth of the closed- loop system resulting in a faster response, while it also 
implements integrating action providing the steady-state accuracy. Although there is some 
conservatism introduced by the fact that upper and lower bounds are used, this example 
demonstrates that a significant increase of the performance can be obtained using this 
method. 



Remark 5.1 In this example the ability of the method to minimize the contribution of the 
slow mode using an additional objective function p{z) in Problem 2.1 was demonstrated. 
Because of the specific objective function (52) the slow mode was completely cancelled 
within the controller. Therefore, for this specific example another way to arrive at con- 
troller (55) is to specify only two closed-loop poles ps^4^ = —2 ± Aj (without pi^2 = ^2j ), 
resulting in C{s) = i^^±^^ which indeed is the minimal form of ([55]). Interestingly, the 



optimization problem results in this controller in an automated and systematic manner. 
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Figure 5: Results exponential bounds relaxation, (a) New (solid black) and original (solid 
grey) step responses, original bounds (dashed grey), and new bounds (dashed black), 
(b) Bode diagrams of the original (grey) and the new (black) controller. 



5.2 Using the multivariate polynomial relaxation 

In this section, the multivariate polynomial relaxation from Section [3l2 is used to obtain 
suitable values of qo,qi,q2 to improve the step response of closed-loop system (43). We 
have that m = 1, = 0, rzc = 2, ai = —1, 02 = —2, /3i = 2 and P2 = 4. Let 6 = 1 such that 
M = cos(r) and f = sin(r) so that (27) yields 

y(t) = ((ai + jbi) {u + jvf + (ai - jhi) {u + jvf) X 



+ {{(12 + Jb2) {u + jv f + (a2 - J62) {u + jv)^) 

= (2ai [u'^-v^] +2bi2uv) X 

+ (2a2 {u'^ + v^-Gu^v^) + 2bi (W- W)) 

As a consequence, for this example we obtain 

Wi{u, v) = — f ^, Tiiu, v) = 2uv 

W2{u, v) = u'^ + v"^ — Gu^v'^, r2{u, v) = Avu^ — Auv^. 

yielding the time response 

y{u, V, X) =yo + {2ai(u^ — t>^) -|- 4:biuv)X 

+ (2a2(M^ + v^- Gu'^v^) + 8b2ivu^ - uv'^))X'^ 



(5^ 



(59) 



(60) 



which is a multivariate polynomial with 3 independent variables {u, v, X) and three decision 
variables (go, Qi, (l2)- Note that (m, v, X) G J-'originai- Since J-'origmai is not the finite union of a 
basic semialgebraic set, we can use Algorithm[T]or use the precomputed overapproximation 
in Section 3.2 to obtain an e-close overapproximation J-'approx of J-'originai- We use here the 
precomputed overapproximation J-'approx with e = e~^'^'" ~ 0.009. In accordance with 
Section |4| we formulate the problem as to find q such that the overshoot 7 is small and 
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that the steady-state error of the step response is minimized. Therefore, the problem is 
posed as 

min 10(1 — yo)"^ + 7 

diel) (61) 



'J0,'?l,iJ2 

s.t. 



- y{u,V,X) >0 V(u,t;, A) G J'approx- 

Rewriting this optimization problem gives 

10(l-yo)' + 7 
(|46| 

'J - y{u,v,\) >0 y{u,v,\)eJ'o 

7 — y{u, f , A) > V(u, V, A) G J^i 

J - y{u,v,X) >0 W{u,v,\) e J^2, 



mm 

go ,91, 92 

S.t. 



(62) 



with J-Q i as in (35) and J-2 as in (33). This optimization problem is then solved with 
a hierarchy of LMI relaxations, as explained at the end of Section 3.2 The size of the 
resulting LMI problem depends on the order of the relaxation, and this can be used as a 
tuning knob to adjust the trade-off between the desired accuracy and the computational 
complexity, see Section [2j Although it is a priori not clear what is the order of the 
relaxation to arrive at the global minimum of 7, a heuristic method is to increase the 
order until not much improvement in the relaxed optimum 7 is observed anymore or 
until one is satisfied with the obtained value of 7. The obtained minimum values of 7 
for various orders of relaxation are given in Table [TJ We used this heuristic approach 
for illustration purposes only. It is recommended to use the more systematic approach 
that is implemented in GloptiPoly 18 to arrive at the global minimum, and certify it 



numerically. Based on the figures in this table we expect that the global minimum of 7 
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297.170 


1.235 


1.235 


1.0718 


1.0718 


1.0718 



Table 1: Upper bound 7 for various orders of LMI relaxation. 

is equal to 1.0718, representing an overshoot of 7.18%. The corresponding Youla-Kucera 
parameter is given by 



q{s) 



which yields the controller 



C{s) 



-32.0 - 17.0607s - 3.0227s' 



3.0s3 + 20.0s2 + 49s + 100 
s3 + 2.0s2 + 10.9s ■ 



(63) 



(64) 



The step responses of both the original closed loop with the c?- minimal controller ( 40 ) and 



of the closed loop with controller (64) are depicted in Fig. |6| which shows a significant 
improvement as expected. The maximum of the step response y{t) equals 1.0714 (i.e., 
7.14% overshoot), which indeed is e-close to 7 = 1.0718. 

This example showed that after controller design by pole placement it is possible to shape 
the transient time response of the system by assigning zeros to the closed-loop system 
through a suitable extension of the controller under pole invariance. 
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1 2 3 4 5 

Time [s] 

Figure 6: New (black) and original (grey) step responses. 

6 Conclusions 

In this paper we provided a generally applicable design framework to obtain linear con- 
trollers for linear systems subject to time-domain constraints. In order to arrive at this 
design framework we extended recent results in [19j that applied only in case of real 
closed-loop poles and external signals (e.g. references or disturbances) having Laplace 
transforms with real poles. The design method is based on synthesizing linear controller 
via a closed-loop pole placement method in which the additional design freedom in terms 
of the Youla-Kucera parameter is used to satisfy time-domain constraints on the closed- 
loop signals based on sum-of-square techniques. In order to extend the method to the 
practically relevant case where both inputs and closed- loop systems with complex conju- 
gate poles are allowed, we proposed two relaxations. 

The first relaxation, called the exponential bounds relaxation, exploited exponential up- 
per and lower bounds on the response to any Laplace transformable input. Although 
this gives rise to potential conservatism, an example showed that by prescribing poly- 
nomial time-domain bounds, the system's performance to a step input can be improved 
with respect to the settling-time and the steady-state error. The second relaxation, the 
multivariate polynomial relaxation, removed the potential conservatism completely as we 
formally proved that it can approximate the original problem with arbitrary accuracy. 
Using these relaxations, we indicated how a polynomial optimization problem could be 
set up, in which next to constraint satisfaction, we could also optimize certain important 
closed-loop properties such as overshoot, settling-time and steady-state error. The result- 
ing optimization problem can be solved using sum-of-squares and convex programming 
methods. As a consequence, the provided design framework is systematic in nature, as 
was also illustrated using a numerical example. 
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